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Abstract. We have re-examined the classical problem of the macroscopic equation of state for a hydrostatic 
isothermal self-gravitating gas cloud bounded by an external medium at constant pressure. We have obtained 
analytical conditions for its equilibrium and stability without imposing any specific shape and symmetry to the 
cloud density distribution. The equilibrium condition can be stated in the form of an upper limit to the cloud 
mass; this is found to be inversely proportional to the power 3/2 of a form factor n characterizing the shape of the 
cloud. In this respect, the spherical solution, associated with the maximum value of the form factor, fi = 1, turns 
out to correspond to the shape that is most difficult to realize. Surprisingly, the condition that defines the onset 
of the Bonnor instability (or gravothermal catastrophe) can be cast in the form of an upper limit to the density 
contrast within the cloud that is independent of the cloud shape. We have then carried out a similar analysis in 
the two-dimensional case of infinite cylinders, without assuming axisymmetry. The results obtained in this paper 
generalize well-known results available for spherical or axisymmetric cylindrical isothermal clouds that have had 
wide astrophysical applications, especially in the study of the interstellar medium. 
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1. Introduction 

Consider a static isothermal gas cloud embedded in a 
medium characterized by an external pressure p^. Under 
equilibrium conditions, it is expected to satisfy the virial 
theorem: 

2K + W = 3 Ph V , (1) 

where K is the total internal kinetic energy associated 
with the thermal motions of the gas molecules and W 
is the total gravitational energy. From the relation K = 
3Mcf/2, we see that 

VhV = Mcl + \w<^. (2) 
6 m 

Here c s is the thermal speed inside the cloud, given by 
c s = \J kT/m, where k is the Boltzmann constant, T the 
cloud temperature, and m the mean mass of the gas par- 
ticles. Hence, the presence of gravity softens the standard 
equation of state pbV = MkT/m for a gas cloud of volume 
V (Boyle's law) into one characterized by an effective tem- 
perature T e g < T. Clearly, T e ff depends on the size, the 
shape, and the temperature of the cloud. In particular, 
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for given values of all other parameters, we could imag- 
ine density distributions for which W is a large negative 
quantity, so that T c g < 0; these distributions would not 
be compatible with the assumed equilibrium conditions. 
On the other hand, even when T Q g > 0, the equilibrium 
configuration might not be viable because it could be un- 
stable. 

Within spherical symmetry, this line of thought has led 
to the remarkable discovery (Bonnor 1956; see also Ebert 
1955) that a thermodynamic type of instability takes place 
when the density contrast (between the center and the 
periphery of the cloud) exceeds a certain threshold value. 
This relatively simple conclusion has opened the way to a 
variety of interesting "applications" , best exemplified by 
the interpretation of the structure of weakly collisional 
stellar systems (globular clusters) in terms of the onset of 
the so-called gravitational catastrophe (see Lynden-Bell & 
Wood 1968 and following papers; in particular, see Spitzer 
1987). The underlying reason for the instability can be 
traced to a basic property of self-gravity, of being asso- 
ciated with a negative specific heat. It would be concep- 
tually interesting to study whether and how the general 
conclusions change for different cloud shapes. 

A study of such well-posed mathematical problem ap- 
pears to be timely. Indeed, there has been enormous 
progress in the study of the cloudy structure of the inter- 
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stellar medium (e.g., Falgarone et al. 1992; Falgarone et al. 
1998; Lada et al. 1994; Alves et al. 1998). Curiously, some 
of the properties that are often considered as analytically 
useful in idealized theoretical models have been shown to 
be, at least in some cases, also realistic. In particular, while 
the possibility of non-trivial shapes for gas clouds is clearly 
suggested by observations of the cold interstellar medium 
(see, e.g., Loren 1989a, Lada et al. 1991), in some cases 
these observations quantitatively support a simple inter- 
pretation in terms of pure hydrostatic equilibrium (see, 
e.g., Alves et al. 2001). Similarly, it has long been rec- 
ognized that cold clouds generally conform to the virial 
equilibrium conditions (e.g., see Myers 1983; Loren 1989b; 
Jijina et al. 1999). Furthermore, although magnetic fields 
and rotation are generally accepted to play a significant 
role in the dynamics of gas clouds (see Goodman et al. 
1995), there is evidence that cloud shapes are decoupled 
from the presence of magnetic fields (compare the case of 
the Taurus cloud, Heyer et al. 1987, with that of the p 
Ophiuchi cloud, Goodman et al. 1990) and that rotation 
can often be insignificant (e.g., Walker et al. 1993); more- 
over, statistical arguments applied to the distribution of 
measured axial ratios for cloud cores suggest that those 
structures are intrinsically prolate (David & Verschueren 
1987; Myers et al. 1991). 

Therefore, we find it natural to address the issue of the 
stability of an isothermal cloud in hydrostatic equilibrium, 
bounded by an external medium characterized by a spec- 
ified constant pressure p^, in its most general case, that 
is beyond the well-known framework of spherical symme- 
try. Since we do not consider the role of magnetic fields 
and, especially, non-thermal motions, we limit our discus- 
sion to the scale of cloud cores, typically a few tenths of 
a parsec. Still, a study on such scales is very important if 
we wish to understand the process of star formation. In 
contrast with previous analytical investigations, we then 
let the shape and density distribution of the cloud devi- 
ate significantly from the spherical shape. The origin of 
these shapes should be traced to the presence of external 
forces (tidal fields and/or other boundary forces). Here we 
will not be concerned about how a specific shape can be 
realized, appealing to the empirical fact that observations 
demonstrate that cloud cores find the way to settle in a va- 
riety of shapes. Note that filamentary structures are very 
common (e.g., Loren 1989a; Bally et al. 1987), and often 
each filament contains several distinct subcondensations, 
sometimes periodically spaced (Dutrey et al. 1991). 

By generalizing the study by Bonnor (1956), in this 
paper we obtain an expression for the compressibility 
dpb/dV at constant temperature and fixed shape and 
thus reconstruct the isothermal curves in the (V,pb) plane. 
These are parameterized by a dimensionless form factor 
p, a quantity that depends only on the shape of the cloud, 
but otherwise follow the qualitative behavior characteris- 
tic of the self-gravitating isothermal sphere; in particu- 
lar, no equilibrium is available above a certain value of 
the boundary pressure. Our analysis quantifies how such 
equilibrium condition depends on p, and shows that the 



sphere, being associated with the maximum value for the 
form factor (/i = 1), corresponds to the shape for which 
equilibrium is most difficult to realize. We then proceed to 
address the condition for the Bonnor instability and show 
that it can be cast in a form that is actually independent 
of the cloud shape. 

The paper is organized as follows. The relevant basic 
relations are introduced in Sect. |^. In Sect. |^ we derive the 
equivalent of Boyle's law for self-gravitating clouds. The 
implications of the results obtained for the equilibrium 
and stability of the cloud are clarified in Sect. |4| and are 
summarized in two simple conditions (Eqs. ( |34|) and (j35|)). 
In Sect. [| we present analogous results for (infinite) cylin- 
ders. Finally, in Sect. |(| we briefly summarize the main 
results obtained in this paper. 

2. Basic relations 

We summarize here the basic relations that will be used 
through this paper. 

Let us call p(x) the mass density distribution and p(x) 
the (scalar) pressure of a cloud at point x. If the cloud is 
in hydrostatic equilibrium, then 

X7p = -p\7<P , (3) 

where 3>(x) is the Newtonian gravitational potential, de- 
scribed, inside the cloud, by the Poisson equation 

V 2 $ = 4nGp . (4) 

Furthermore, if the cloud is isothermal we can write 

p = c\p . (5) 

Taking T = 10 K and a cloud made of molecular hydrogen, 
we find a typical value of c s ~ 200 m s _1 . 

Equations (||-|5|) combined together lead to Emden's 
equation (Emden 1907): 



(6) 



where po is a characteristic density (e.g., the maximum 
density of the cloud). A useful property of Eq. (g) that is 
the basis of the following analysis is the homology theorem 
(Chandrasekhar 1967): If p(x) is a solution of this equa- 
tion, then X 2 p(Xx) is also a solution. In fact, we can write 
Eq. (JT|) in dimensionless form by introducing the variables 

«=£ (7) 
Po 

and £ = sx, where s^ 1 is a suitable "Jeans length" (cf. 
Jeans 1929) 



UnGpo 



(8) 



(9) 



Hence u is solution of the equation 

V|ln«(0 = -u(e), 

where V| is the Laplacian with respect to the variable £ 
In the following we will reserve the subscript b for quan- 
tities measured at the boundary of the cloud. 
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3. Boyle's law 

In this section we will derive Boyle's law for a self- 
gravitating cloud of arbitrary shape subject to a uniform 
external pressure pb- In particular, we will obtain an ex- 
pression for the derivative dpb/dV at constant tempera- 
ture and fixed shape, for a cloud of given mass. A similar 
calculation had been performed by Bonnor (1956; see also 
Ebert 1955) in the case of a spherical cloud. Here we imag- 
ine that the cloud is embedded in an external medium in 
an arbitrary dynamical state. Our main assumption is that 
such an arbitrary state enforces a non-trivial geometry for 
the cloud, while the cloud has had time to settle into a hy- 
drostatic equilibrium and is thus described by Eq. (||) ■ The 
possibility of such non-trivial shapes is clearly suggested 
by observations of the cold interstellar medium (see, e.g., 
Loren 1989a, Lada et al. 1991), which in some cases also 
support the hydrostatic hypothesis (see, e.g., Alves et al. 
2001). 

Calculations will be carried out using the following 
scheme: We will first calculate the change of the cloud 
volume corresponding to a change of the density param- 
eter po', then we will relate the variation of the density 
parameter po to the variation of the pressure pb at the 
boundary of the cloud. The shape of the cloud is taken to 
remain unchanged. 

Since the cloud is in equilibrium with an external pres- 
sure pb, the boundary of the cloud is defined by an iso- 
density surface characterized by pb = Pb/c^. When the 
external pressure changes, the cloud boundary remains an 
iso-density surface. In other words, the transformations 
considered involve only a change of po (and thus of s) in 
p(x), while remains unchanged. Finally, we note that, 
since the boundary is an iso-density surface, Vp and Vjit 
are normal to the cloud boundary. 

3.1. Perturbations at constant mass 

Suppose that the density parameter po changes by a quan- 
tity Spo. Then the cloud density distribution will change 

by 

ds 

8p(x) = Spou(sx) + 6p - — poVfutsx) ■ x . (10) 
dpo 

As a result, the mass inside a given volume V will change 

by 



SM{V) = / 5p{x)dV 
Jv 



(11) 



From Eq. (|l^) we obtain 

5M(V) = Sp [ u(sx) dV 
Jv 

ds f 

+ Spo- — po / Veu(sx) ■ x dV . 
dpo Jv 



(12) 



The first term on the r.h.s. of the previous equation is 
simply SpoM(V)/po. For the second term, we can apply 



the general relation V(/t>) = V/ • v + /V • v. Thus we 
obtain 



SM(V) 



S Po M(V) 
Po 

— u(sx)V ■ x] dV 

Sp M(V) r dins 

\-Opo-; 

Po dp 



3 / p u(sx)dV 



. ds po f r— / , \ \ 

6po- / [V • (u(sx)x) 

upo s Jv 



d> pou(sx)x ■ n dS 

JdV 



(13) 



where we have used the Gauss theorem (dV is the bound- 
ary of V, n is the unit vector associated with the oriented 
surface element). Since p is constant on dV, we find 



SM(V) = 



5p Q M(V) | 
Po 

Sp M{V) 



8po 
2po 



Pb 



[ V -xdV -3M(V) 
Jv 



35p p h V 



(14) 



2p 2p 

Therefore, for perturbations at constant cloud mass a vari- 
ation Spo of the "central" density must be compensated 
for by a suitable change of volume 5V . Since the change 
of mass due to a change of volume SV is pb<5^, we must 
have 

5p M 3(5p pbV 



5Vp h = o 



(15) 



2p 2p 

This equation controls the change of central density Spo 
when the volume of the cloud changes. 

3.2. Total mass in virial equilibrium 

For the following calculations it is useful to obtain a sim- 
ple alternative expression for the total mass applicable to 
a self-gravitating isothermal cloud in hydrostatic equilib- 
rium. For the purpose, we will use Eq. (^) in the integral 
of the mass: 



dV 



4nG 



Vp • n 



dS . 



dv 



(16) 



This expression can be simplified by noting that p is con- 
stant on dV and Vp • n = — ||Vp||, because the boundary 
is defined as an iso-density surface. [Here ||A|| = V A ■ A 
denotes the norm of the vector A] Hence, we find 



AnGpb J dv < 



V ( u\\ dS . 



(17) 



The use of this relation is equivalent to the use of the 
equation for virial equilibrium. 

3.3. Compressibility 

We now evaluate the change of the external pressure of 
a cloud as the cloud volume changes, i.e. the deriva- 
tive dpb/dV at constant temperature, constant mass, and 
fixed shape. 
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If we allow for changes of size, at the boundary Eq. ( |10| 
becomes 



5pt 



Po 



ds 

po- — dpo\/£U(sx) ■ x 
0-Po 

p sSV 



\V iU (sx)\\ dS 



(18) 



TdV 



The last term on the right-hand-side of this equation can 
be explained as follows. Since the cloud can change its 
size, we expect a related variation of the external density 
Pb- Such change can be written as 5p = — V/o(x) \\5x for 
each point x, where 6x is the "stretch" of the boundary 
at x (we recall that this stretch is always parallel to Vp) . 
We also note that Sp is the same for all points on the 
boundary of the cloud, because by definition this quantity 
can be identified with Spb in the case where Spo = 0. Thus 
the volume changes by 



dS 



,)T = </> SxdS = -Sp I ^ 
hv Jdv 



(19) 



Inverting this equation, we obtain that the change of den- 
sity due to a modification of the size of the cloud is given 

by 



5~ P = - 



SV 



p sSV 



iVpW^dS 



7av 



§dv\\ V t u ( sx )\\ dS 



(20) 



which is the last term of Eq. (|18|) . 

The second term on the r.h.s. of Eq. (|l8|) is constant 
for any point x on the boundary dV because all the other 
terms in the equation are so. Thus we can average this con- 
stant term on the boundary using a convenient weighted 
mean: 



V^u(sx) ■ x 



V^u(sx) ■ x 
av \\V^u(sx)\\ 

x ■ ndS J 



dS 



dS 



av 



ov ||Vfu(sx) 
dS 



ov 



\V^u(sx)\ 



§ av \\X7 ( u(sx)\\ dS 



(21) 



With this device the second term is reduced to a form that 
is similar to that of the last term of Eq. (|l8|). 

3.4. Generalized Boyle's law 

Let us briefly summarize the main results obtained so far. 
We wish to obtain an expression for dpb/dV at constant 
temperature and mass or, equivalently, an expression for 
Spb as a function of SV. Equation ( p"8| ) is close to our 
need, but, unfortunately, it contains the auxiliary quantity 
5pQ. On the other hand, we can use Eq. ( |l5| ) to express 
Spo in terms of SV and thus eliminate Spo from Eq. (|l8|). 
Following this procedure, we obtain the desired dpb/dV, 
but the resulting expression contains integrals that are, at 
first sight, difficult to interpret. 

We now note that, with the help of Eq. (|l7|), a key 
term involving an integration of ||Vfu|| (see Eqs. ( |l8[ ) 



and can actually be recognized to be proportional 

to GM/y 4 ' 3 : 



[V f w(sx) dS 



TdV 



4tt\ 1/3 MGppb 



(22) 



This expression is reminiscent of what enters in the dis- 
cussion of the spherical cloud by Bonnor (1956). For the 
purpose, we have introduced the dimensionless form factor 
as follows: 



12tt 



3 \ 1 / 3 F 4 /3 



47T 



s 2 \s 



dV 



|V e u(sx)|| dS 



— (p || V^u(sa;)|| 1 dS 
& Jdv 



(23) 



In this way, we have isolated the contribution from inte- 
grals of ||V^u||. We postpone a description of the form 
factor to the next subsection, where the reason for the nu- 
merical factors adopted in its definition will become ap- 
parent. 

Using relation ( p2[ ) in Eq. ( ^l| ) and combining it with 
Eqs. (f|) and |[|) we find 



5pb = - 



3p h V - M 



, 47T 

1 -«T 



1/3 



GpM 2 
6c 2 p h V 4 / 3 



(24) 



Finally, using Eq. (w we obtain the desired result: 



dp b 
dV 



2 Ph 



4tt\ 1/3 GpM 2 



T.M 



1 



(25) 



This expression generalizes Eq. (2.16) in the article by 
Bonnor (1956). Before discussing the physical interpreta- 
tion of the equation just found and its consequences, the 
form factor p deserves a special digression. 

3.5. Form factor 

The form factor, which plays a key role in the equilibrium 
and stability of the cloud, is defined in terms of integrals 
of ||V^m|| and ||V^u|| on the boundary dV . Actually, 
the complex definition ( p3| ) has a simple interpretation in 
terms of general properties of p. 

3.5.1. p: is dimensionless 

As already noted, p is a dimensionless quantity, which is a 
choice that allows us to preserve the scaling of the terms 
as found in Bonnor's (1956) derivation. In this respect, 
note that the term V 4 ^ 3 compensates for the two surface 
integrals in Eq. (|23|). As shown below, the introduction 
of a dimensionless quantity allows us to characterize and 
summarize the properties associated only with the shape 
of the cloud. 
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3.5.2. n is scale invariant 

It is easily verified that the form factor fi is scale invariant, 
i.e. it does not change under the transformation x i— > tx, 
with t > 0. In other words, the form factor depends on 
the shape of the cloud but does not depend on its size. 
Actually, the scale invariance is directly related to the fact 
that fj, is dimensionless and to the vectorial notation used 
in its definition. 



3.5.3. /j, = 1 for a sphere 

In the case of a sphere, the form factor reduces to unity. 
In fact, in this case by symmetry ||V^u(sa;)|| is constant 
on the surface of the sphere, so that 



\i = 12tt 



4tt 



1/3 



^4/3 



1 . 



(26) 



Here S is the area of the surface of the sphere. Thus, from 



Eq. (25), we recover Bonnor's (1956) result for the sphere. 
Note that the numerical factors in the definition of \i have 
been chosen precisely to recover fi = 1 for the sphere. 




4 6 
V (10 46 m 3 ) 

Fig. 1. The p-V plot for an isothermal spherical cloud 
made of molecular hydrogen at T = 10 K (c s ~ 
200 m s^ 1 ) with mass M = 1 Mq. The isothermal lies 
on a simple pV — const curve at large volumes, but then 
at smaller sizes departs from it and the cloud becomes un- 
stable. The region of stability is on the right side of point 
A (which corresponds to the maximum of the boundary 
pressure pb)- The solution spirals into the critical point C. 



3.5.4. n never exceeds unity 

The form factor can be thought of as the product of two 
terms, /x = ^1^2, with (cf. Eq. (pq)) 



Ml = 12n [ - 



3 \ 1/3 yi/3 



S 2 



(27) 



The first term, fj,i, compares the volume V of the cloud 
with its surface S. This quantity equals unity when the 
cloud is a sphere and is always smaller than unity in the 
non-spherical case. In fact, the sphere is the solid which 
has the largest volume at given surface (this intuitive the- 
orem has a far from trivial proof, due to De Giorgi 1958). 
The second term, \i2 , is the ratio between the generalized 
harmonic mean of Vgit on the boundary dV and the 
generalized simple mean of the same quantity. The har- 
monic mean is always smaller than the simple mean, the 
two being equal only when all the values involved are con- 
stant. Thus fj,2 < 1 and fi% = 1 if and only if V^u is 
constant on dV (this happens for the isothermal sphere). 
In conclusion, /i < 1 in general and [L = 1 only for the 
sphere. 

4. The Bonnor instability 



Equation ( |25| ) allows us to discuss the Bonnor instability 
(often known as gravothermal catastrophe; see Lynden- 
Bcll & Wood 1968) for clouds of arbitrary shape. This 
instability, we recall, is due to a change in the sign of the 
derivative of pb with respect to V. If this derivative is neg- 
ative, as normally happens when gravity can be neglected, 
a slight decrease of the volume will produce an increase 
of the cloud boundary pressure, which would tend to re- 
store the initial configuration. If, instead, this derivative is 



positive, then a small decrease in the volume of the cloud 
would correspond to a reduction of the internal boundary 
pressure and thus to a collapse of the cloud. 

In order to illustrate the instability process in more 
detail, let us first refer to the spherical case. We have 
integrated the differential equation ( |25| ) numerically under 
the boundary condition pbV — c 2 M for V — > 00. Figure [l] 
shows the case of a spherical cloud of molecular hydrogen 
at T = 10 K (c s ~ 200 m s" 1 ) with mass M = 1 M Q . The 
plot clearly shows that the cloud follows approximately 
the standard law pV = const at large volumes, but when 
sufficiently compressed it exhibits the Bonnor instability. 
In particular, the p-V curve has a maximum and then 
spirals into a singular point. All points of the p-V curve 
where the derivative is positive correspond to unstable 
equilibria, since a small decrease of the volume will make 
the cloud collapse. Bonnor has shown that actually, in the 
case of a spherical cloud, all points in the p-V plot between 
A and the critical point correspond to unstable equilibria 
(below we will generalize this result to all clouds). The 
idea is that a spherical cloud with pressure and volume 
beyond A will present an instability for compressions of 
some internal "core" (we will discuss this instability in the 
general case below). 

In the general case of an arbitrary shape, the situation 
turns out to be similar. In fact, going back to Eq. (p5|), 
we can see that the form factor enters the global equation 
of state together with the gravitational constant G. In 
other words, a non-spherical cloud behaves like a spher- 
ical cloud with a reduced gra vitational constant Gfi. As 
discussed above (see Sect. 3J5), the isothermal sphere has 
the largest admissible value for the form factor, fi = 1, and 
thus for the spherical shape the effects introduced by self- 
gravity are largest. The Gfi invariance also implies that the 
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Fig. 2. The p-V plot (in logarithmic coordinates) for 
isothermal clouds of different shapes made of molecular 
hydrogen at T = 10 K with mass M = 1 M . The form 
factors (U are, from top left to bottom right, 0.125, 0.25, 
0.5, and 1. In logarithmic coordinates, a change of p is 
described by a simple translation of the isothermal curve. 



isothermal p-V curves exhibit similar shapes regardless of 
the value of p. In Fig. |we have plotted the integral curves 
of Eq. ( p5| ) under the assumption that the cloud form fac- 
tor is constant as V changes. Figure || clearly shows that 
the Bonnor instability occurs at any p > 0; moreover, 
the effect of reducing the form factor is to postpone the 
onset of the instability to smaller volumes. In the limit 
of a vanishing form factor we recover the standard law 
pV = const. From this figure we also see that in loga- 
rithmic coordinates a change of the form factor p simply 
results in a simple translation of the p-V curve. In fact, it 
is easily checked that, if we define 

P = pp\ V = Vp-\ (28) 

then Eq. (|2^) is transformed back into the equation ap- 
propriate for a sphere, since p drops out. As a result, a 
curve relative to a form factor p can be obtained from the 
curve for the sphere by multiplying the volume by p 3 and 
the pressure by p~ 3 . 

Point A is a zero of dp^/dV and thus satisfies 

{AvV'^GpM 2 , . 

PA=[ — ) -±IJT ■ ( 29 ) 



or 



Using the Emden (1907) solution for the sphere, and tak- 
ing advantage of the p 3 scaling, we can give the position 
of point A for arbitrary shapes: 

3 



V A ^ 0.290307 



PA 



1.39766 



G 3 p 3 M 2 



(30) 



(31) 



The curves in the p-V plane spiral into the critical point 
C. This point is easily obtained by imposing that both the 
numerator and the denominator of Eq. ( p5| ) vanish. Then 
we obtain 



7T / GpMY 



6 V 



PC = — 



7T G 3 p 3 M 2 



(32) 
(33) 



Let us now complete the discussion of the Bonnor in- 
stability for a general cloud. We have already noted that 
all clouds for which the derivative dp^/dV is negative are 
unstable. Actually, we can show that all points in the p-V 
curve between A and C are unstable. In fact, suppose that 
the cloud is in a state represented by a point B between 
A and C, and let us consider smaller and smaller portions 
of the cloud with boundaries characterized by a constant 
pressure. Let us call V(p) and M(p), respectively, the vol- 
ume and the mass of the part of the cloud with pressure 
larger than p. We now observe that a small compression of 
this part of the cloud is still described by Eq. J2q) with the 
obvious replacements pb l— * P, V i— > V(p), and M i— > M(p). 
As we increase the value of p, we approach the central part 
of the cloud. Assuming that the maximum cloud density 
po is finite, we have M(V) < poV, and the numerator 
of the last term of Eq. fl25|) will tend to unity because 
M 2/y±/3 < ^2^2/3 ^ o as U -> 0. As a result, in the 
central part of the cloud the gas follows the standard law 
pV = const, i.e. we are in the right part of the plot p-V. If 
we now move back from the central part to the boundary 
of the cloud represented by point B we must travel along 
the p-V curve (actually, we are moving along different p- 
V curves, with respect to the reference one, because the 
mass M(V) of the sub-cloud changes while we consider 
larger and larger volumes) and we will reach, at a certain 
boundary pressure p a , a point of type A. In other words, 
we have shown that the sub-cloud characterized by p > p a 
is on a point of type A of the p-V plot. If this sub-cloud 
has a small negative fluctuation of its volume V a = V(p a ), 
a decrease of its boundary pressure will occur and this 
will lead to a further reduction of its volume. Thus we 
must conclude that every cloud characterized by a point 
in the p-V plot between A and C is unstable. This is the 
generalization of the Bonnor instability to non-spherical 
clouds. 

Point A then defines the boundary of the stability re- 
gion. If the external pressure is larger than pa, no equilib- 
rium solution exists for the cloud. This provides an upper 
limit for the cloud mass: 



M < M r 



1.18223- 



y/ G 3 p 3 p cxt ' 



(34) 



where p cxt is the external pressure (that we suppose to be 
balanced by the cloud boundary pressure ph)- Inserting 
here a typical value for the pressure of the interstellar 
medium, p cxt = 1.4 x 10 -13 Pa, and using c s = 200 m s _1 
we find M crit ^ 4.7p~ 3/2 M Q . 
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If condition (|4|) is met, we can still be in the instability 
region if we are in the spiral part of the p-V curve. This 
happens if pbV < paYa, and thus a necessary condition 
for stability, independent of p, is 



P < Pc 



2.46457 



Pcxt 



(35) 



where p — M/V is the mean density of the cloud. 
Inserting here the typical values used before, we find 
p cr it ~ 8.6 x 10~ 18 kg m~ 3 , corresponding to about 2600 
molecules per cm 3 . Equation ([55|) can also be read in terms 
of a density contrast requirement p < 2.46457pt,. Equation 
( |35| ) proves that the condition for the Bonnor instability 
can be cast in a form that is independent of the shape of 
the cloud. However, Bonnor's stability result for spheres, 
when expressed as a condition po/pb < 14.0, does not hold 
for clouds of arbitrary shape. 

Finally, we wish to emphasize that Eqs. (HJ) and ( ^5| ) 
are only necessary conditions for equilibrium and stability, 
and not sufficient conditions. In other words, other kinds 
of instabilities may occur for clouds that satisfy Eqs. ( |34| ) 
and ([Bo]). Moreover, the presence of rotation or magnetic 
fields may qualitatively change these results, as suggested 
by recent investigations on the equilibrium of clouds in 
the presence of magnetic fields (see Fiege & Pudritz 2000; 
Curry & Stahler 2001). 

5. Boyle's law for (infinite) cylinders 

So far we have assumed that the cloud has finite volume. 
However, several observed clouds appear as filaments and 
can be well modeled by cylinders. The technique described 
above, suitably adapted, can be used also to obtain Boyle's 
law for (infinite) cylindrical clouds of arbitrary shape, de- 
scribed by a mass density of the form 



p(xi,x 2 ,x 3 ) = p(x 1 ,x 2 ) 



(36) 



Note that no assumption on the function p(xi,x 2 ) has 
been made, and in particular we do not assume a circular 
section. We can now calculate the compressibility dpb / dV 
following the steps taken earlier in this paper. Since the 
calculations are very similar to the 3D case, we will just 
outline the main differences here. 

Equations (§) to (0) remain unchanged; however, 
since all functions are independent of x 3l it is conve- 
nient to define here x — (xi,x 2 ), £ = (^1,^2)1 and 
V = (d/dxi,d/dx 2 ). The "total" mass cannot be defined 
for infinite cylinders, and hence we will refer to a linear 
mass density A defined as 



A= / p(x)dS 
Is 



(37) 



where S is the section of the cylinder. 



5.1. Perturbations at constant mass 

Using Eq. ©, which remains unchanged in the 2D case, 
we can write the change of "mass" SX(S) when the central 



density changes from po to po + Spo assuming no variation 
in the cloud section S. The result obtained is basically the 



same as in Sect. 3.1 with the important difference that now 
V • x — 2. As a result, Eq. (|l4|) now becomes 



6X(S) = 



SpaX(S) 5p 



Po 
SpoPbS 

Po 



+ 



2po 



Ph 



V • x dS — 2\(S) 



(38) 



Note that the two terms involving X(S) cancel each other. 
If we now allow for a change of section 5S and impose 
that the total "mass" A be constant we find the analogue 
ofEq. @: 



Spo 
Po 



s 



(39) 



This equation is independent of pb and simply implies 
po oc 1/5. 

5.2. Compressibility 

Equation (|l8|) rewritten for cylinders is 



5ph = 



Spopb 
Po 



ds 

Po j — OpoV^SX) ■ X 

dpo 
PosSS 



\V^u(sx)\\ dt 



(40) 



JdS 



Similarly to the 3D case, we can rewrite the second term 
on the r.h.s. in a more convenient way by using a suitable 
average. In this case we have 



V^u(sx) ■ x 



V^u(sx) ■ x 
OS \\V£u(sx)\\ 
2S 



df 



df 



as ||Vfu(sx)| 



ros 



|V s u(sx) d£ 



(41) 



Inserting this result in Eq. (|4(j) and using relation (|3^) we 
find 



5pb = - 



SSpb 
S 



(42) 



This simple result should be compared with Eq. (|24|) , valid 
in the three-dimensional case. We finally find the analogue 
OfEq. (H), 



dpb 
dS 



Pb 

' s 



(43) 



This equation states that a self- gravitating isothermal 
cylinder behaves like a normal gas under compressions, 
i.e. it follows the standard law p^S = const. However, 
we should emphasize that in general the constant that ap- 
pears on the r.h.s. of this expression is not equal to Ac;?. 
In other words, an isothermal cylinder at temperature T 
follows a law similar to the one of a perfect gas but at 
an ef fect ive temperature T c ff < T (see Eq. (^); see also 
Sect. 5^ below). The difference T — T c g will depend on 
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the particular shape of the cylinder and on its linear mass 
density A. Furthermore, for cylinders, since no form factor 
needs to be introduced, Eq. ([l3|) can be integrated exactly 
without additional assumptions. 

Since the cylinder follows a law similar to the standard 
law for perfect gases, no Bonnor instability is present in 
this case. Still we have to make sure that T e g > 0. From 
Eq. ( ff3| ) we see that this is true as long as an equilibrium 
exists with p^ > 0. In the following subsection we will 
consider an axisymmetric cylinder in order to illustrate 
this point by means of a simple example. 

5.3. Axisymmetric cylinder 

The case of a circular cylinder, with density of the form 
p(x) = p(r) dependent only on the distance r from the 
axis, has already been the subject of several investiga- 
tions (see, e.g., Ostriker 1964; Horedt 1986; Bastien et 
al. 1991). For an axisymmetric cylinder the exact form of 
the solution of Emdcn's equation is known (Stodolkiewicz 
1963; Ostriker 1964). In our notation, the solution can be 
written as 



1. The generalized Boyle's law for an isothermal self- 
gravitating cloud in hydrostatic equilibrium has been 
studied. No assumption has been made on the shape 
of the cloud. 

2. An expression for the compressibility dpt,/dV has been 
obtained and, from that, two conditions for stable equi- 
librium have been derived. 

3. The equilibrium of the cloud has been shown to be con- 
trolled by one dimensionless form factor. The spherical 
shape turns out to be the one for which the effect of 
self-gravity is largest. 

4. A simple relationship between the isothermal curves in 
the p-V plane and the form factor has been obtained. 

5. The actual condition for the Bonnor instability can be 
written in a form that does not depend on the shape 
of the cloud. 

6. For infinite cylinders of arbitrary shape we have shown 
that the Bonnor instability is suppressed. 

Acknowledgements. We would like to thank Joao Alves for 
helpful discussions. Part of this work has been carried out at 
the Scuola Normale Superiore, Pisa, Italy. 



p(r) = Pq 



1 + 



(sr) 2 



(44) 



We can integrate this expression out to a certain radius 
R in order to get the linear density A of the cylinder. On 
the other hand, if we know the linear density A and the 
radius R of the cloud, we can obtain the central density 
Pq. A simple calculation yields 



A 



Po 



S 



GXS 
~2cT 



(45) 



As anticipated in Eq. (|39|), we find po oc 1/S. Inserting 
this in Eq. (^J) at r — R and using Eq. (||) we finally find 



p h S = Ac s 1 



GX 

2c? 



(46) 



This is a standard isothermal law, equivalent to Boyle's 
law, but characterized by an "effective" temperature 



T cS = T 



GXr 



2k 



(47) 



The term GXm/2k represents the contribution from the 
potential energy of the cylinder, as indicated by Eq. (||). 
From Eq. ( |46| ) we immediately see that when the cloud lin- 
ear density is larger than a critical density A cr it = 2c 2 s /G, 
the cloud cannot be in equilibrium (the required boundary 
pressure pb would be negative). 

6. Conclusions 

The main results obtained in this paper can be summa- 
rized as follows. 
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